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The  following  are  summaries  of  various  projects  completed  for  the  Air  Force  project: 

1.  Intercalation  of  Carboxymethyl-P-Cyclodextrin  into  Magnesium-aluminum  Layered 
Double  Hydroxide  (published  in  J.  Chem.  Soc.  Dalton  Trans.) 

New  intercalated  compounds  have  been  prepared  by  incorporating  carboxymethyl-|3- 
cyclodextrin  (CMCD)  into  magnesium  aluminum  layered  double  hydroxide  (Mg/Al  LDH).  The 
sorptive  uptake  of  CMCD(3)  and  CMCD(14),  with  a  degree  of  carboxymethyl  substitution  of  j 
and  14,  respectively,  by  Mg/Al  LDH  was  examined  at  65  °c,  and  the  resulting  complexes  were 
characterized  using  XRD  and  FT-IR.  Results  indicated  that  Mg/Al  LDH  could  uptake 
approximately  twice  as  much  CMCD(3)  as  compared  to  CMCD(14).  Results  of  XRD  and  FT-IR 
confirmed  that  both  CMCD(14)  and  CMCD(3)  could  be  intercalated  into  Mg/Al  LDH  interlayers 
due  to  an  ion  exchange  process.  After  exposure  to  CMCD(14)  and  CMCD(j),  the  d-spacing  of 
Mg/Al  LDH  expanded  from  8.74  A  to  15.48  A  and  20.63  A,  respectively.  Intercalated 
CMCD(14)  molecules  formed  monolayer  coverage  in  the  Mg/Al  LDH  interlayers,  with  the 
cavity  axis  perpendicular  to  the  LDH  layer,  while  CMCD(3)  molecules  adopted  either  a  parallel 
monolayered  arrangement  or  a  perpendicular  bilayer-like  coverage.  Both  CMCD(14)  and 
CMCD(3)  molecules  were  believed  to  be  loosely  packed  within  the  Mg/Al  LDH  interlayers. 

2.  Molecular  Inclusion  Characteristics  of  Hydrophobic  Organic  Compounds  by  a 
Modified  B-cyclodextrin  Inercalated  Within  a  Layered  Double  Hydroxide  (Revised 
and  resubmitted,  should  be  accepted  by  J.  Inch  Phenom.  Mol.  Recog.  Chem.) 

A  study  was  conducted  to  evaluate  the  inclusion  properties  of  various  nonionic 
hydrophobic  organic  compounds  by  a  novel  intercalate  derived  from  magnesium-aluminum 
layered  double  hydroxide  (Mg/Al  LDH)  and  carboxymethyl-p-cyclodextrin  with  a  degree  of 
substitution  of  3  [CMCD(3)j.  The  isotherm  sorption  results  at  25  °c  showed  that  the  CMCD(3)- 
Mg/Al  LDH  intercalate  could  retain  all  the  organic  compounds  studied  and  its  sorption  affinity 
for  organic  compounds  was  positively  related  to  their  hydrophobicities.  The  host-guest 
interaction  was  attributed  to  a  partition  process  of  the  organic  compounds  into  the  cyclodextrin 
cavity,  while  a  stereoselective  interaction  might  also  be  involved  due  to  the  intercalation  of 
CMCD(3)  within  Mg/Al  LDH  interlayers. 
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3.  Selectivity  and  Molecular  Sieving  Effects  of  Organic  Compounds  on  the  Intercalate 
of  B-cycIodextrin  and  Layered  Double  Hydroxide  (Submitted  to  Clays  and  Clay 
Minerals) 

This  study  dealt  with  the  sorption  characteristics  of  10  nonionic  hydrophobic  organic 
compounds  by  a  novel  intercalate  derived  from  magnesium-aluminum  layered  double  hydroxide 
(Mg/Al  LDH)  and  carboxymethyl-P-cyclodextrin  with  a  degree  of  substitution  of  14 
[CMCD(14)].  The  isotherm  sorption  results  at  25  °C  showed  that  the  CMCD(14)-Mg/Al  LDH 
intercalate  exhibited  size/shape  selectivity  for  the  organic  compounds  studied.  This  selectivity 
appeared  to  result  from  the  intermolecular  micropores  and  the  hindered  accessibility  of  the 
cyclodextrin  cavity.  Such  selectivity  and  sieving  effect  may  be  useful  in  chromatographic 
applications  as  well  as  in  selective  removal  of  organic  pollutants  from  water. 

4.  BTEX  Sorption  by  Organo-Clays :  Cosorptive  Enhancement  and  Equivalence  of 
Interlayer  Complexes  (Published  in  Soil  Sci.  Soc  Am.  J.) 

Aqueous  organic  contaminants  are  effectively  sorbed  by  organo-clays.  Organo-clays 
were  prepared  from  montmorillonite  (SWy-1,  SAz-1)  and  vermiculite  (VSC)  clay  minerals  by 
exchanging  quaternary,  alkylammonium  cations  for  the  inorganic  exchange  cations,  thus,  a 
variety  of  clay-charge/organic-cation  combinations  were  formed.  Dodecyl-trimethylammonium 
(DTMA),  hexadecyltrimethylammonium  (HDTMA),  didodecyl-dimethylammonium  (DDDMA), 
and  dioctadecyldimethylammonium  (DODMA)  organic  cations  were  used.  An  aqueous  mixture 
of  BTEX  (benzene,  toluene,  ethylbenzene,  xylenes)  consistent  with  unleaded  gasoline  was  used 
to  relate  organo-clay  sorptivity  to  structure.  Greater  sorption  of  the  BTEX  mixture  constituents 
than  the  pure  compounds  occurred  when  total  sorbed  BTEX  loadings  on  the  organo-clays 
exceeded  10  g/kg.  A  co-sorption  effect  caused  curvilinear  isotherms  and  enhanced  BTEX 
mixture  sorption;  BTEX  sorption  increased  the  organo-clay  organic-matter  content  which 
promoted  additional  BTEX  sorption. 

The  sorptivities  of  the  organo-clays  were  generally  proportional  to  the  organic  carbon 
contents;  although  SAz-DDDMA  and  VSC-HDTMA,  with  20  to  30%  less  carbon,  had  twice  the 
BTEX  sorptivity  of  SAz-DODMA  and  VSC-DDDMA.  No  relationship  was  found  between 
surface  area  and  BTEX  sorption.  A  straight-chain  alkyl  cation  conferred  5  times  greater  BTEX 
sorption  to  an  organo-clay  than  did  a  comparable  cyclic  alkyl  cation.  This  suggests  that 
rearrangements  in  the  alkyl  groups  may  facilitate  sorption  in  a  manner  analogous  to  solvation. 
Organo-clay  BTEX  sorptivity  was  also  governed  by  the  type  of  alkylammonium  complex;  a 
paraffin  complex  with  a  25.4  to  27. 6A  basal  spacing  was  optimal.  Organo-clays  of  similar  basal 
spacing  prepared  using  different  clays  and  alkylammonium  cations  had  equivalent  sorptive 
phases.  Selection  of  an  appropriate  combination  of  clay  and  organic  cation  can  optimize 
contaminant  sorption  by  organo-clays. 
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5.  Sorption  of  the  Ionizable  Organic  Compound,  Dicamba,  (3,6-Dichloro-O-Anisic  Acid) 

by  Organo-Clays  (Published  in  Chemosphere) 

The  sorption/desorption  of  3,6-dichloro-o-anisic  acid  (dicamba)  by  hexadecyl-trimethyl- 
ammonium  (HDTMA)-  and  dioctadecyldimethylammonium  (DODMA)-exchanged  smectites 
was  studied  at  different  concentrations  and  pH  levels  to  evaluate  the  capability  of  organo-clays  to 
sorb  an  ionizable  hydrophobic  compound.  Dicamba-organo-clay  complexes  were  also  prepared 
and  characterized  by  x-ray  diffraction  (XRD)  and  Fourier  transform  infrared  spectroscopy 
(FTIR)  to  determine  the  possible  sorption  mechanisms  involved.  The  results  indicated  that 
dicamba,  although  water-soluble,  ionizable,  and  potentially  mobile,  was  strongly  sorbed  by  the 
organo-clays  and  was  not  readily  desorbed.  Sorption  isotherms  could  be  divided  into  two 
distinct  regions;  a  non-linear  region  at  low  dicamba  concentrations  and  a  linear  region  at  high 
dicamba  concentrations;  partitioning  was  believed  to  be  the  dominant  process  involved  in  the 
sorption  of  dicamba  at  higher  concentrations. 

Organo-clays  with  higher  organic  carbon  contents  (i.e.,  amount  and  molecular  weight  of 
exchanged  organic  cation  (HDTMA  or  DODMA))  usually  resulted  in  greater  dicamba  sorption. 
Solution  pH  significantly  affected  the  sorption  of  dicamba  by  organo-clays,  with  the  isotherm 
inflection  point  near  the  dicamba  pKa  value.  Nearly  twice  as  much  dicamba  could  be  sorbed  in 
the  molecular  form  as  compared  to  its  anionic,  deprotonated  form.  Excessive  DODMA  was 
found  to  increase  dicamba  sorption  at  pH  levels  greater  than  the  dicamba  pKa,  which  was 
possibly  caused  by  charge  reversal  of  the  organo-clay  surface.  Possible  sorption  mechanisms, 
based  on  XRD  and  FT-IR  data,  suggest  hydrophobic  and  ionic  attractions  both  contributed  to  the 
sorption  of  dicamba. 

6.  Sorption  of  Trichloroethylene  by  Organo-clays  in  the  Presence  of  Background 

Humic  Substances  (Submitted  to  Water  Research) 

A  study  examined  trichloroethylene  sorption  characteristics  by  organo-clays 
(hexadecyltrimethylammonium  (HDTMA)  and  didodecyldimethylammonium  (DDDMA) 
exchanged  smectites  (SWy-1  and  SAz-1))  in  aqueous  solutions  both  in  the  presence  and  absence 
of  humic  substances  (aldrich  humic  acid  and  peat  fiilvic  acid),  as  well  as  the  impacts  of  the 
preloading  of  humic  substances.  The  results  indicated  uptake  of  humic  substances  by  organo- 
clays  and  the  effects  of  humic  substances  on  trichloroethylene  sorption  by  organo-clays  was 
dependent  on  the  type  of  organo-clay  and  humic  substance.  Selection  of  an  appropriate 
combination  of  clay  and  organic  cation  can  optimize  the  removal  of  trichloroethylene  as  well  as 
humic  substances  from  aqueous  solutions  by  organo-clays.  Organo-clays  may  be  an  alternative 
to  activated  carbon  sorbents  in  specific  water  treatment  processes. 
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7.  Removal  of  Polynuclear  Aromatic  Hydrocarbons  by  Orgairo-clays  from 

Water/methanol  Mixtures  (In  preparation  for  Environ  Sci.  and  Technol.) 

A  study  was  conducted  to  evaluate  the  effectiveness  of  10  organo-clays  for  removal  of 
PAHs  (pyrene  and  anthracene)  in  the  presence  of  methanol.  Preliminary  results  indicated  that 
cosolvent  theory  could  be  successfully  used  to  predict  the  sorption  behaviors  of  PAHs  by  organo- 
clays  from  water/methanol  mixtures,  and  that  organo-clays  had  strong  affinities  for  both  pyrene 
and  anthracene.  However,  the  sorption  affinities  varied  among  the  organo-clays,  depending  on 
the  combination  of  clay  and  organic  cation.  This  study  has  implications  on  understanding  the 
usefulness  of  organo-clays  as  sorbents  for  PAHs  pollutants  in  real-world,  complex  waste  streams 
such  as  industrial  wastes  and  landfill  leachates  which  are  usually  composed  of  a  mixture  of 
water,  pollutants,  and  miscible  organic  solvents. 

8.  Modeling  BTEX Sorption  by  Organo-Clay  Derivatives  (In  preparation  for  Soil  Sci.  Soc 
Am.  J.) 

The  objectives  of  the  present  study  were  to  quantify  BTEX  (benzene,  toluene, 
ethylbenzene,  o-,  m-,  p-  xylene)  sorption  by  various  organo-clay  derivatives  using  different 
isotherm  models  and  to  evaluate  competitive  sorption  characteristics  of  these  hydrocarbons  in 
relation  to  organic-cation  mass  and  organo-clay  basal  spacing.  Two  montmorillonite-derivatives 
including  trimethylphenylammonium  (WNM-C9H14N+)  and  trimethylammonium  adamantane 
( WNM-C 1 3H24N+),  and  three  hectorite-derivatives  including  trimethylphenylammonium  (CNH- 
C9HI4N+),  trimethyl-ammonium  biphenyl  (CNH-C15H18N+)  and  trimrthylammonium  fluorene 
(CNH-C,6HlgN+)  were  used  in  this  study.  The  organo-clays  were  prepared  by  exchanging  the 
inorganic  cations  in  the  clay  minerals  with  the  quaternary  alkyl  ammonium  cations.  A  batch 
sorption  study  was  conducted  with  a  BTEX  mixture  that  contained  the  constituent  compounds  in 
proportions  consistent  with  the  composition  of  unleaded  gasoline. 

The  results  for  all  of  the  six  compounds  indicated  that,  the  concentrations  in  sorbed  phase 
increased  curvilinearly  with  that  in  solution  phase.  Corresponding  to  the  equilibrium 
concentrations  the  amounts  of  sorbed  hydrocarbons  were  higher  for  C9H|4N  derivatives  as 
compared  to  CI3H24N+  derivative  followed  by  CI5H18N  andC|6H18N  derivatives.  Solid  phase 
partitioning  of  ethylbenzene,  o-xylene  and  m-xylene  was  maximum  in  CNH-C9H14N\  and  that  of 
benzene  was  greatest  in  WNM-C9H14N+.  Hydrocarbon  sorption  data  were  statistically  (p  <  0.1) 
fitted  to  Langmuir,  Freundlich  and  Temkin  models.  The  Langmuir  model  was  fitted  to  BTEX 
sorption  on  WNM-C9H14N+,  WNM-C13H24N+  and  CNH-C9H14N*  whereas  adsorption  of  CNH- 
C15HI8N+  and  CNH-C16H18N+  did  not  conform  to  this  model.  The  site-densities  of  WNM- 
C9H14N+,  WNM-C  13H24N+  and  CNH-C9H14N+  for  BTEX  were  estimated  to  be  1.89-1 1.4x1 017/m2, 
2.80-19.1xl017/m2,  and  1. 60-44.2x1 017/m2,  respectively,  indicating  a  high  order  (1017)  of 
magnitude  for  all  three  organo-clays.  Relatively  higher  values  of  the  adsorbate  buffering 
capacity  (K,B,)  of  the  C9H14N+  derivatives  than  WNM-C13H24N+  suggested  that  sorption  on 
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C9HI4N+-clays  would  be  more  effective.  The  K,B,  further  indicated  that  the  values  for  the 
montmorillonite  derivative  were  either  comparable  to  or  higher  than  the  hectorite  derivative  for 
the  same  cation  C9HI4N+.  The  binding  affinity  (Kf)  values  obtained  from  the  Freundlich  model 
was  consistently  higher  for  C9H|4N+  derivatives  as  compared  to  the  WNM-Ci3H24N  .  However, 
the  maximum  Kf  values  were  noted  for  CNH-Ci5H18N+  and  followed  by  CNH-C16H18N+.  The 
ranges  of  nf  values  for  benzene,  toluene,  ethylbenzene,  o-xylene,  m-xylene  and  p-xylene  were 
0.90-1.75,  0.78-1.74,  0.64-1.88,  0.61-1.62,  0.65-1.67,  and  0.63-1.88,  respectively.  The  lowest  nf 
values  for  all  six  hydrocarbons  were  consistently  observed  for  CNH-Ci5HigN  ,  and  the  highest 
values  were  mostly  found  in  case  of  C9H|4N+-clays.  The  Temkin  constant  Ktl,  pertaining  to  heat 
of  chemisorption,  was  found  to  decrease  linearly  with  molecular  weight  of  the  cations  present  in 
the  organo-clay  derivatives  thereby  suggesting  the  organo-clays  with  low  molecular  weight 
cations  would  have  greater  heat  of  chemisorption,  which  physically  means  greater  hydrocarbon 
retention. 

The  slope  of  a  differential  change  in  K„  with  cationic  molecular  weight  followed  the 
order  of  toluene  >  m-xylene  >  p-xylene  >  o-xylene  >  ethylbenzene  >  benzene.  Greater  Kom 
values  indicated  greater  partitioning  of  the  hydrocarbons  into  the  organo-clays;  thus,  maximum 
partitioning  of  all  six  hydrocarbons  will  be  caused  by  WNM-C9H14N+  organo-clay  followed  by 
CNH-C9Hi4N+  and  WNM-C13H24N+  derivatives;  CNH-C15HI8N+  and  CNH-C16H18N+  organo-clays 
corresponded  to  much  lower  Kom  values;  a  similar  sequence  in  the  magnitude  of  aG  (sorption 
free  energy)  also  supported  this  finding.  The  positive  (+)  values  of  log  Kom/Kow  ratio  indicated 
that  WNM-C9HI4N+,  WNM-C13H24N+  and  CNH-C9H,4N+  derivatives  could  result  in  greater 
partitioning  of  BTEX  as  compared  to  octanol  medium.  Preferential  partitioning  of  benzene  and 
toluene  was  also  indicated  by  the  greater  magnitude  of  log  Kom/Kow,  which  was  due  to  greater 
solubility  of  these  two  hydrocarbons  as  compared  to  ethylbenzene  and  the  xylenes.  This  was 
further  reflected  in  the  plot  showing  the  effect  of  organic-cation  thickness  on  BTEX  partitioning. 
Corresponding  to  the  d001  values,  log  Kom/Kow  for  benzene  and  toluene  were  higher  than  other 
hydrocarbons.  The  log  Kom/Kow  values  for  all  six  hydrocarbons  decreased  following  the  same 
pattern  as  the  basal  spacing  increased.  This  means  greater  partitioning  of  the  hydrocarbons  will 
occur  with  organo-clays  with  lower  interlayer  spacing,  for  example  C9HI4N+  derivatives. 
Knowledge  obtained  from  this  study  is  expected  to  contribute  to  the  understanding  of  future 
contaminant  remediation  techniques. 
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may  serve  as  a  substrate  for  adsorption  of  non-polar  species.  Associated  with  the  adsorption 
phenomenon,  organic  anions  accelerate  the  dissolution  of  these  silicates,  so  the  presence  of  such 
organic  compounds  may  affect  the  physical  and  chemical  properties  of  the  aquifer. 
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The  primary  focus  of  our  component  of  the  effort  has  been  the  development  of  in  situ 
spectroscopic  methods  to  study  the  adsorption  and  desorption  of  various  contaminant  model 
compounds  at  mineral  surfaces.  Toward  this  end,  we  have  developed  a  new  methodology  for 
using  ATR-FTIR  (attenuated  total  reflection-Fourier  transform  infrared)  spectroscopy  to 
probe  adsorption  of  various  organic  compounds  at  mineral  surfaces.  In  this  method,  a  flow 
cell  is  constructed  so  that  it  can  be  placed  against  an  ATR  optical  element  (such  as  a  ZnSe 
crystal)  in  such  a  way  as  to  trap  the  mineral  particles  within  close  proximity  to  the  surface  of 
the  ATR  element.  Due  to  the  evanescent  optical  field  that  exists  just  outside  of  the  ATR 
element,  the  infrared  spectrum  of  the  portion  of  the  particles  that  are  within  the  evanescent 
region  can  be  sampled,  as  well  as  the  spectra  of  any  adsorbed  species  that  are  present.  In  this 
way,  one  obtains  the  vibrational  spectrum  of  the  adsorbate  species. 

There  are  several  intrinsic  advantages  to  this  technique.  The  most  significant  one  is 
that  the  strong  interactions  that  lead  to  adsorption  of  most  adsorbates  will  produce  significant 
changes  in  the  vibrational  spectrum  of  the  adsorbate,  allowing  it  to  be  distinguished  spectrally 
from  the  corresponding  species  in  solution.  Thus,  this  method  will  typically  provide  the 
capability  to  differentiate  easily  between  solution  phase  species  and  adsorbed  species. 
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Figure  1  -  Schematic  of  the  ATR-FTIR  flow  cell  used  to  study  adsorption  at  the  surfaces  of  mineral 

particles. 

This  method  has  been  applied  to  several  different  systems.  Our  initial  experiments 
employed  silica  particles  as  the  adsorbent  phase.  In  one  set  of  experiments,  we  examined  the 
adsorption  of  species  bearing  carboxylate  functionalities  (e.g.  acetate)  at  silica  surfaces.  As 
can  be  seen  below  in  Figure  2,  we  found  that  the  adsorption  of  such  species  does  not  occur 
spontaneously.  This  is  because  over  the  pH  range  that  we  studied  (4-8),  the  surface  charge 
of  silica  is  neutral  to  slightly  negative,  leading  to  electrostatic  repulsion  of  the  negatively 
charged  acetate.  Further,  the  silica  surface  is  not  expected  to  provide  sufficiently  electrophilic 
adsorption  sites  for  acetate  (i.e.  there  are  no  attractive  sites  at  which  acetate  could  adsorb  via 
complex  formation).  However,  when  metal  ions,  such  as  Ca2\  are  present  in  solution  which 
can  enter  into  complexes  with  the  surface  siloxane  groups,  the  carboxylate  species  adsorb 
onto  the  surface.  This  adsorption  is  presumably  due  to  complexation  between  the  adsorbed 
metal  ions  and  the  carboxylates  functionalities.  A  manuscript  on  this  study  is  in  preparation. 


Upper:  Acetate  adsorbed  on  Si02  in  presence  of  Ca2+ 


Figure  2  -  Spectra  of  acetate  (Ac)  anion  in  the  pore  volume  of  SiO:  particles.  Upper:  0. 1  M  NaAc  + 
lffiuM  CaCb  .  Lower:  0. 1  M  NaAc  without  CaCk  Obtained  in  ATR  mode  with  ZnSe  element  on  a 

Bomem  MB  100  FTIR 


The  lower  curve  in  the  Figure  2  shows  the  spectrum  of  acetate  anion  in  aqueous 
solution  (peaks  at  1560  cm'1  from  the  antisymmetric  carboxylate  stretch  and  1415  cm'1  from 
the  symmetric  carboxylate  stretch),  these  bands  are  observed  because  of  the  portion  of  the 
acetate  solution  that  is  trapped  within  the  pore  volume  of  the  Si02  particles  adjacent  to  the 
ZnSe  ATR  element.  The  upper  curve  shows  the  spectrum  that  is  obtained  when  the  SiO: 
particles  are  exposed  simultaneously  to  acetate  and  Ca  ion.  This  spectrum  clearly  shows  the 
responses  from  both  the  solution-phase  acetate  and  the  adsorbed  acetate,  as  indicated  in  the 
figure.  The  peaks  from  the  adsorbed  acetate  (shoulder  at  1540  cm'1  and  peak  at  1495  cm'1) 
are  shifted  from  those  of  the  solution-phase  species  in  a  manner  that  is  consistent  with 
bidentate  interaction  with  the  adsorbed  Ca2'  (i.e.  both  acetate  oxygens  are  attached  to  the 
adsorbed  acetate). 

These  data  were  important  because  they  validated  the  essential  assumptions  behind 
the  choice  of  the  ATR-FT1R  technique,  namely  that  it  would  be  possible  to  spectrally 
differentiate  the  solution  phase  and  adsorbed  species  by  virtue  of  changes  in  their  vibrational 
spectra.  A  manuscript  on  this  study  is  in  preparation. 

The  attached  manuscript  (submitted  to  Langmuir)  describes  a  study  of  the  adsorption 
of  the  surfactant  dodecyltrimethylammonium  (DTAB)  cation  at  silica  surfaces.  These  systems 
are  of  importance  to  the  (JSAF  because  they  have  the  potential  to  be  used  in  remediation  of 
subsurface  contaminants,  such  as  TCE.  Figure  3  below  shows  the  isotherm  obtained  for 
adsorption  of  DTAB  at  the  silica  surface,  along  with  isotherms  calculated  for  various 


adsorption  models.  These  data  demonstrate  convincingly  that  it  is  also  possible  to  observe 
adsorbed  species  even  when  their  vibrational  spectra  are  not  substantially  different  from  those 
of  the  solution  phase  species. 


DTAB  isotherm  on  Si02 


•  Experimental  data 

-  Frumkin  isotherm  fit 

Langmuir  isotherm  fit 
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Figure  3  -  Adsorption  isotherm  for  DTAB  at  silica,  with  calculated  model  isotherms,  demonstrating 

Frumkin-like  adsorption  behavior. 

We  have  achieved  several  long  term  goals  that  derive  from  this  general  area  of 
interest.  First,  we  have  developed  a  method  that  can  be  used  to  identify  the  nature  of  the 
interactions  between  the  contaminants  and  the  mineral  surfaces.  Our  objective  in  this  regard 
was  to  develop  an  ability  to  predict  the  extent  of  adsorption  of  a  given  contaminant  on  a 
given  mineral,  where  the  prediction  will  be  based  solely  on  the  molecular  structure  of  the 
contaminant  and  the  nature  of  the  mineral  surface.  Second,  we  developed  the  experimental 
methodology  that  will  allow  us  to  measure  the  desorption  rates  of  contaminants  from  mineral 
surfaces.  The  objective  of  this  component  of  the  work  was  to  provide  a  set  of  kinetic  data  on 
desorption  rates  that  will  be  of  use  to  that  part  of  the  scientific  community  involved  in 
modeling  the  subsurface  transport  of  contaminants.  The  reasoning  behind  this  effort  is  to  gam 
a  better  understanding  of  dispersive  effects  in  subsurface  transport.  Third,  we  investigated  the 
adsorption  of  amphiphilic  and  hydrophobic  compounds.  The  objective  here  was  to  learn  how 
the  adsorption  of  surfactants  influences  the  adsorption  and  desorption  of  hydrophobic 
contaminants. 


An  In-situ  ATR-FTIR  Study  of  DTAB  Adsorption  at  a  Silica  Gel  Surface 

Hong  Gu  and  Daniel  A.  Buttry* 

Department  of  Chemistry 
University  of  Wyoming 
Laramie,  WY  82071-3838 
BUTTRY@UWYO.EDU 


Introduction 

The  adsorption  of  surfactants  from  solution  at  the  solid-liquid  interface  is  of  great 
importance  both  because  of  theoretical  interest  in  such  processes  and  because  of  the  many 
applications  in  which  they  are  important.  As  a  consequence  of  this  high  level  of  interest,  there 
are  many  recent  studies  of  the  phenomenon.1"  Applications  in  which  surfactant  adsorption  is 
important  include  detergency,  ore  floatation,  microelectronics  and  environmental 
remediation.  A  particularly  interesting  area  has  to  do  with  the  use  of  surfactants  in 
environmental  remediation.4'5  There  are  at  least  two  distinct  modes  of  operation  in  such 
efforts.  In  one,  the  micellar  phase  of  the  surfactant  is  used  to  mobilize  hydrophobic 
subsurface  contaminants,  so  that  they  may  be  removed  more  effectively  using  pump-and-treat 
methods.  In  this  application,  surfactant  adsorption  is  not  desirable,  because  it  reduces  the 
amount  of  surfactant  that  is  available  to  form  micelles.  In  another  mode  of  operation,  the 
surfactant  is  used  to  temporarily  immobilize  a  hydrophobic  subsurface  contaminant  by 
trapping  it  within  the  hydrophobic  layers  created  by  surfactant  adsorption  at  mineral  surfaces. 
In  this  case,  adsorption  of  the  surfactant  is  critical  to  the  success  of  the  temporary 
immobilization  Thus,  information  about  the  adsorption  of  surfactants  at  mineral  surfaces  is 
essential  for  the  proper  engineering  of  such  remediation  efforts.  Due  to  the  usually  negatively 
charged  nature  of  sediments  and  soil  particles,  cationic  surfactants  such  as  DTAB  can  adsorb 
strongly  onto  these  particles.  However,  little  is  known  about  the  details  of  such  adsorption  at 
the  microscopic  level.  This  study  describes  the  use  of  in-situ  ATR-FTIR  methods  to  monitor 
the  adsorption  of  DTAB  at  the  aqueous  solution/silica  interface. 

The  use  of  in-situ  FTIR  to  study  the  solid/liquid  interface  has  a  history  that  can  be 
traced  back  to  the  early  1970’s.  Rochester  and  co-workers6'7  studied  adsorption  of  organic^ 
molecules  from  carbon  tetrachloride  solutions  onto  silica.  Tejedor-Tejedor  and  co-workers 
applied  the  cylindrical  internal  reflection  (CTR-FTIR)  method  to  the  study  of  adsorption  at 
the  mineral  oxide/water  interface.  By  utilring  the  attenuated  total  reflection  (ATR)  method, 
some  mineral/water  interface  have  been  studied,  such  as  the  silicaAvater  interface  , 
goethite/water  interface8'9,  alumina/water  interface1'12,  etc.  These  studies  showed  that  ATR  is 
capable  of  semi-quantitative  characterization  of  adsorption  behavior  at  aqueous  interfaces.  A 
principal  advantage  of  this  method  over  other  methods  is  that  the  limited  penetration  depth  of 
the  beam  into  the  solution  reduces  the  absorption  of  the  beam  by  water,  which  is  the  usual 
problem  with  infrared  experiments  in  aqueous  solutions.  Comparison  to  ex-situ  FTIR 
methods  used  in  surfaces  studies'^’14  shows  that  another  advantage  of  the  m-situ  ATR-FTIR 


method  is  that  it  has  the  possibility  to  be  used  to  monitor  the  kinetics  of  the  adsorption  and/or 
desorption  processes. 

Experimental  Materials  and  Methods 

Materials.  Dodecyltrimethylammonium  bromide  (CnFEs^CFE^^Br )  was  purchased  from 
Aldrich  with  purity  greater  than  99%  and  used  as  received.  Silica  gel  for  column 
chromatography  was  obtained  from  Aldrich  and  used  without  further  purification.  This  silica 
gel  has  a  BET  surface  area  of  500  m2/g  and  pore  volume  of  0.75  cm3/g.  The  distilled  H20 
used  for  the  aqueous  solutions  was  deionized  and  filtered  by  using  a  0.45  /im  nylon 
membrane. 

FTDR  Spectroscopic  Measurements.  Attenuated  total  reflection  IR  spectra  were  taken 
using  a  BOMEM  Michelson  MB- 100  FTER  spectrometer  with  a  liquid  nitrogen  cooled 
HgCdTe  detector.  ZnSe  internal  reflection  elements  were  from  Harrick  Scientific  Co.  and 

had  dimensions  of  50x1 0x3 mm  and  employed  a  45  angle  of  incidence.  The  optical  bench  of 
the  spectrometer  was  purged  with  dry  nitrogen.  Spectra  of  the  samples  were  obtained  with 
256  scans  at  a  resolution  of  2  cm  1  Figure  1  shows  the  flow  cell  developed  and  used  in  our 
laboratory  to  flow  solutions  past  the  silica  particles.  This  flow  cell  is  designed  to  fit  into  a 
standard  ATR  sample  holder.  Typically,  the  silica  gel  particles  are  introduced  into  the  body  of 
the  flow  cell  and  contained  within  the  sample  compartment  (which  comprises  the  space 
between  the  membrane  and  the  ZnSe  waveguide)  by  a  membrane  fabricated  from  a  thin  sheet 
of  A1  into  which  ca.  50  pm  holes  had  been  drilled.  This  membrane  constrains  the  particles  to 
press  aaainst  the  ZnSe  waveguide.  Then,  the  flow  cell  is  pressed  against  a  silicone  gasket, 
and  therefore,  against  the  ZnSe  waveguide  using  a  screw  press.  In  this  configuration,  the 
silica  particles  are  pressed  against  the  waveguide,  and  therefore,  are  within  the  evanescent 
wave  region  of  the  IR  beam. 

Before  any  spectra  were  taken,  the  cell  with  dry  silica  was  placed  into  the  sample 
chamber,  which  was  then  purged  with  dry  nitrogen  for  10  minutes.  Then,  a  2.7mM  NaCl 
solution  was  pumped  into  the  cell  for  10  minutes  at  a  flow  rate  of  2.5  mL/min.,  and  the  first 
single-beam  spectrum  was  taken.  This  spectrum  was  used  as  the  reference  for  all  the  sample 
spectra  thereafter.  Preliminary  experiments  had  shown  that  an  equilibrium  could  be  reached 
at  the  silica  gel  surface  after  1 0  minutes  of  exposure  to  a  DT AB  solution  within  the 
concentration  range  of  0.5  mM  ~  20  mM.  Thus,  each  sample  spectrum  was  taken  after 
pumping  a  DTAB  solution  of  a  known  concentration  through  the  flow  cell  for  10  minutes  at 
a  flow  rate  of  2.5  mL/min.  After  the  spectrum  had  been  taken,  a  2.7  mM  NaCl  solution  was 
pumped  into  the  cell  for  10  minutes  to  wash  out  the  adsorbed  DTAB.  Note  that  even  the 
DTAB  solutions  contained  2.7  mM  NaCl.  The  purpose  of  the  use  of  low  concentration  salt 
solutions  in  these  experiments  was  to  provide  for  the  charge  compensation  that  is  necessary 
whenever  a  charged  species  such  as  DTA*  adsorbs  or  desorbs. 


Results  and  Discussion 

Figure  2a  shows  the  spectrum  of  a  1  mM  aqueous  solution  of  DTAB  at  a  bare  ZnSe 
ATR  element.  This  spectrum  shows  the  CH  stretching  envelope  that  is  characteristic  of 
species  with  long  alkyl  chains.  It  also  shows  an  envelope  in  the  1500  -  1700  cm'1  range.  This 


envelope  is  due  to  water  absorption  from  the  OH  bending  vibration,  with  small  absorptions 
from  water  vapor  in  the  spectrometer  chamber  superimposed  on  it.  The  appearance  of  this 
envelope  is  due  to  small  changes  in  the  amount  of  water  within  the  region  of  solution  probed 
by  the  beam.  These  changes  arise  both  from  changes  in  the  mole  fraction  of  water  within  the 
probed  region  between  the  background  and  the  sample  scans  and  from  changes  in  the 
penetration  depth  between  these  two  scans.  This  latter  effect  is  caused  by  subtle  changes  in 
the  refractive  index  of  the  solution  due  to  changes  in  the  solution  composition.  A  full 
description  of  these  effects  will  be  reported  in  a  forthcoming  full  paper.  The  important 
message  to  be  garnered  from  the  spectrum  in  Figure  2a  is  that  the  spectral  region  from  1500 
to  1700  cm'1  is  relatively  obscured  by  water  features,  which  makes  it  difficult  to  use  it  for 
quantitative  analysis  of  DTAB  adsorption.  Thus,  the  CH  stretching  features  around  2800- 
3000  cm'1  have  been  used  instead,  with  spectral  assignments  made  according  to  published 
data.  15 


Figure  2b  shows  spectra  of  DTAB  of  different  concentrations  in  the  presence  of  silica 
particles.  The  CH2  symmetric  stretch  (vs)  is  observed  at  2855  cm'1.  The  CH2  asymmetric 
stretch  (v„)  is  observed  at  2926cm'1.  The  bands  are  more  clearly  assigned  in  Figure  2c,  which 
is  for  the  20  mM  case.  It  is  expected  that  the  alkyl  chain  CH  stretching  vibrations  should  be 
essentially  identical  for  the  adsorbed  and  solution  phase  DTAB  species  ,  because  the  weak 
physisorption  expected  for  this  species  is  not  likely  to  significantly  perturb  the  alkyl  chains. 
Thus,  we  have  integrated  the  spectral  region  between  2846  and  2875  cm  1  to  quantitatively 
evaluate  both  species  simultaneously. 

Figure  3  shows  the  integrated  areas  of  the  spectral  region  between  2846  and  2875 
cm’1  for  various  concentrations  of  DTAB  adsorbed  on  silica  gel.  Series  1  shows  the 
experimental  data.  These  integrated  intensities  include  the  contributions  from  both  bulk 
solution  and  adsorbed  DTAB.  A  method  used  by  Simpson  and  Harris  was  applied  to 
subtract  the  bulk  solution  contribution  to  the  total  intensity.  Since  the  ATR  intensity  is 
proportional  to  concentration  according  to  Beer’s  law,  18  a  linear  regression  was  made  to 
model  the  solution  phase  behavior  .  By  fitting  the  high  concentration  points  (2.5mM~20mM) 
to  the  linear  model,  the  slope  b  and  intercept  a  could  be  generated.  The  intercept  was  the 
intensity  of  the  spectrum  when  DTAB  reached  saturated  surface  coverage  ,  and  the  slope  was 
the  proportional  factor  in  Beer’s  law.  Thus  the  bulk  solution  contribution  to  the  intensity 
could  be  generated  by  using 

Ablution  =  b*C  0) 

where  C  is  the  concentration  of  the  solution  at  equilibrium. 

Figure  3,  series  2  gives  the  integrated  intensities  after  point-by-point  subtraction  of 
the  solution  contribution.  These  data  clearly  have  the  form  of  an  adsorption  isotherm,  and 
suggest  that  the  surface  is  saturated  with  DTAB  at  concentrations  of  2. 5  -  3.0  mM.  In  order 
to  obtain  a  quantitative  measure  of  the  adsorption  behavior,  both  the  Langmuir  and  Frumkin 
models  were  used  to  fit  this  experimental  adsorption  isotherm.  For  both  models,  the 


fractional  coverage  9  was  obtained  by  normalizing  the  integrated  intensities,  i.e.  each 
integrated  intensity  was  ratioed  to  the  extrapolated  intercept  a.  For  the  Langmuir  model 

0=  KC/(1+  KC)  (2) 

where  K  is  the  adsorption  equilibrium  constant,  and  C  is  the  equilibrium  concentration  of  the 
adsorbate,  which,  in  an  accurate  definition,  should  be  the  activity  of  the  adsorbate.  However, 
given  that  the  data  in  this  case  are  significantly  below  the  critical  micelle  concentration  for 
DTAB,  we  have  chosen  to  use  the  concentrations.  The  result  of  the  Langmuir  fit  is  shown  in 
Figure  4a.  The  curve  fit  gives  K=1.6  x  103  M*1,  using  a  non-linear  least  squares  analysis. 

Comparison  of  the  experimental  data  points  with  the  Langmuir  fit  in  Figure  4  makes  it 
clear  that  the  Langmuir  isotherm  does  not  adequately  describe  the  adsorption  of  DTAB  at 
silica.  Thus,  one  must  consider  the  deficiencies  of  the  Langmuir  model,  which  is  based  on  the 
following  assumptions:  (1)  only  a  monolayer  is  formed  at  the  surface  ,(2)  the  surface  sites  are 
homogeneous,  in  the  sense  that  the  adsorption  energy  is  constant  across  the  surface;  and  (3) 
no  interactions  occur  between  adsorbates.17’19’20'21  Consideration  of  the  experimental  points 
reveals  that  the  first  assumption  is  reasonable;  the  adsorption  shows  clear  saturation  behavior 
characteristic  of  monolayer  formation.  However,  the  validity  of  the  remaining  two 
assumptions  is  less  obvious  by  simple  inspection  of  the  data.  In  order  to  test  whether  the 
inclusion  of  adsorbate  interactions  in  the  model  would  provide  a  better  fit.  a  Frumkin  model 
was  then  used.  In  the  Frumkin  model, 

0/(  1- 0)  =  KCexpfe0)  (3) 

the  isotherm  is  modified  by  an  exponential  factor,  the  physical  meaning  of  which  is  that  the 
adsorption  energy  is  a  function  of  surface  coverage.  If  the  interaction  parameter,  £>0,  the 
interaction  between  the  adsorbates  is  attractive;  if  g<0,  it  is  repulsive. 

The  solid  line  in  Figure  4  gives  the  result  of  a  non-linear,  least  squares  fit  of  the 
experimental  data  to  a  Frumkin  model,  with  K  =  4.5  x  10'  M  1  and  g  =2.4,  which  suggests 
attractive  interactions  between  adsorbed  DTAB  molecules,  probably  mediated  by 
hydrophobic  interactions  between  the  alkyl  chains.  The  regression  coefficient,^,  for  the 
Frumkin  fit  is  r2=0.93,  while  the  corresponding  value  for  the  Langmuir  fit  is  r2=0.72.  Thus, 
the  Frumkin  fit  is  seen  to  be  much  better  than  the  Langmuir  fit. 

Conclusion 

An  in-sit u  FTIR-ATR  method  has  been  described  which  can  be  used  to  quantitatively 
measure  the  adsorption  of  surfactants  at  the  mineral/water  interface.  For  DTAB,  a 
monolayer  is  formed  on  the  silica  gel  surface  at  concentrations  near  2.5  -  3.0  mM.  The 
adsorption  can  be  modeled  as  a  Frumkin  isotherm.  The  shape  of  the  isotherm  suggests  that 
the  adsorption  proceeds  to  saturation  (i.e.  monolayer  formation)  and  then  stops.  The 
positions  of  the  CH  stretching  peak  maxima  suggest  that  the  adsorbed  DTAB  species  form  a 
relatively  disordered  structure  on  the  silica  surface.  The  positive  interaction  parameter,  g  = 
2.4,  can  be  interpreted  as  indicating  that  the  adsorbed  DTAB  molecules  attract  one  another. 


This  value  of  the  Frumkin  interaction  parameter  is  in  relatively  good  agreement  with  a  value, 
g  =  2.1,  reported  by  Peng  and  Rusling21  for  DTAB  adsorption  onto  a  Au  electrode  surface. 
Thus,  even  though  the  surfaces  are  distinctly  different,  this  comparison  reveals  that  the 
interactions  between  the  alkyl  chains  are  quite  similar  in  these  two  cases. 

A  significant  conclusion  of  this  study  is  that  the  in-situ  ATR-FTIR  method  provides  a 
relatively  facile  route  into  studies  of  the  adsorption  of  surfactants  at  the  mineral-water 
interface.  Thus,  it  is  now  possible  to  make  comparable  measurements  on  a  wide  variety  of 
surfactant  and  mineral  surface  systems.  In  addition,  it  should  be  possible  to  use  these  methods 
to  examine  the  interactions  between  these  surfactant  systems  and  the  pollutants  themselves 
25,  an  effort  that  is  ongoing  in  our  laboratories.  Specifically,  we  are  now  examining  the  uptake 
on  hydrophobic  pollutants,  such  as  TCE,  into  such  surfactant  layers  in  an  effort  to  understand 
how  to  optimize  this  uptake. 
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Figure  1  A  side  view  of  the  special  flowing  cell  developed  and  used  in  this  experiment. 
Solution  can  flow  through  the  silica  particles  at  the  ZnSe  and  silica  interface, 
thus  an  in-situ  spectrum  can  be  obtained. 

Figure  2  (a)Spectrum  of  ImM  DTAB  when  scan  time  is  1024  and  resolution  is  2  cm'1. 
From  this  spectrum  CH2  stretches  can  be  seen  at  2800-2900cm'1.  Those 
sharp  peaks  at  1400-1800  cm'1  are  from  water  vapor  in  the  IR  chamber. 

(b)  CH2  stretches  shown  in  various  concentrations  of  DTAB. 

(C)  Spectrum  of  20mM  DTAB  at  the  region  of  2850-3000cm'1.  This  spectrum 
shows  both  CH2  stretches  and  CH3  asymmetric  stretch. 

Figure  3  Spectroscopic  isotherm  of  DTAB.  Series  1  shows  the  integrated  intensity  at 
2855cm'1  plotted  vs.  concentration  of  DTAB.  Series  2  shows  the  isotherm 
after  contribution  of  the  solution  part  to  the  intensity  has  been  subtracted. 


Figure  4  Langmuir  and  Frumkin  fittings  of  the  isotherm.  Both  observation  and  statistic 
test  show  that  Frumkin  model  fits  better. 
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Work  on  this  aspect  of  the  project  began  officially  in  January,  1995.  The  people  involved 
include  the  following: 

•  Myron  B.  Allen,  Professor  of  Mathematics 

•  Stacy  Howard,  Ph.D.  candidate,  Mathematics 

•  Patrick  O’Leary,  Ph.D.  candidate,  Mathematics. 

The  aim  of  the  mathematical  modeling  part  of  the  project  was  to  characterize  the  effects  of 
random,  spatially  varying  adsorption  on  the  transport  of  organic  and  inorganic  groundwater 
contaminants. 

More  specifically,  we  view  a  soil  column  or  a  one-dimensional  porous  medium  as  an  array 
of  adsorption  sites,  with  adsorption  characteristics  obeying  a  specified  frequency  distribution 
and  having  a  specified  spatial  correlation  structure.  This  model  reflects  the  uncertainties 
associated  with  heterogeneous  geologic  media.  Given  this  stochastic  model  of  the  porous 
medium,  we  investigate  whether  one  can  model  the  ensemble- average  contaminant  transport 
usin0'  a  standard,  advection-diffusion-type  transport  equation  with  ” effective”  (upscaled) 
coefficients  to  capture  the  effects  of  stochasticity. 

The  methodology  used  has  been  a  combination  of  numerical  and  analytical  investigation, 
based,  for  the  time  being,  on  the  simple  case  of  equilibrium  adsorption  according  to  a  linear 
isotherm.  We  developed  a  finite-element  code  for  the  advection-diffusion-retardation  equa¬ 
tion  with  spatially  varying  retardation  and  a  statistical  algorithm  for  assigning  retardation 
factors  to  cells  in  the  grid  from  a  probability  distribution  having  specified  mean,  variance, 
and  spatial  covariance.  This  numerical  scheme  shows  clearly  that,  when  the  adsorption 
characteristics  are  spatially  correlated  over  distances  comparable  to  the  dimensions  of  the 
contaminant  plume  (or  larger),  the  spreading  of  the  average  plume  cannot  be  adequately 
modeled  using  the  standard,  Fickian  model  for  hydrodynamic  dispersion. 

The  analytical  studies  have  pursued  several  tacks.  One  of  the  most  fruitful  has  been  small 
perturbation  theory,  in  which  we  derive  approximate  equations  governing  the  average  plume 
by  assuming  that  the  retardation  field  has  the  form  of  a  constant  value  plus  a  small,  spatially 
varying  fluctuation  having  mean  zero.  This  approach  yields  a  nonlocal,  integrodifferential 
equation  for  the  average  contaminant  plume.  While  this  equation  explains  the  non-Fickian 
behavior  observed  in  the  numerical  studies,  its  nonlocal  form  is  extremely  inconvenient  for 
numerical  approximations,  and  therefore  it  is  unlikely  to  find  widespread  use  in  applications. 

The  investigation  has  not  ended,  even  though  the  funding  period  has.  We  are  now 
investigating  numerically  tractable  methods  to  solve  for  the  average  plume  -  methods  that 
do  not  require  the  solution  of  nonlocal  equations. 

The  work  yielded  the  following  publications  and  invited  lectures: 
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P.  O’Leary,  M.B.  Allen,  and  F.  Furtado,  ’’Groundwater  transport  with  stochastic  retarda¬ 
tion,”  Proceedings,  Third  IMACS  Conference  on  Iterative  Methods,  ed.  by  J.  Wang  et  al., 
Jackson,  WY,  9-12  July,  1997. 

J.  Yang,  R.  Zhang,  J.  Wu,  and  M.B.  Allen,  ’’Stochastic  analysis  of  adsorbing  solute  trans¬ 
port  in  three-dimensional,  heterogeneous,  unsaturated  soils,  Water  Resources  Research  33.8 
(1997),  1947-1956. 

M.B.  Allen  and  F.  Furtado,  ’’Computational  methods  for  porous-media  flows,”  Chapter  6, 
Advances  in  Fluid  Mechanics:  Porous-Media  Flows,  ed.  by  J.P.  du  Plessis,  Computational 
Mechanics  Publications,  Southampton,  U.K.,  1997,  255-302. 

J.  Yang,  R.  Zhang,  J.  Wu,  and  M.B.  Allen,  ’’Stochastic  analysis  of  adsorbing  solute  in 
two-dimensional  unsaturated  soils,”  Water  Resources  Research  32:9  (1996),  2747-2756. 

B.  Liu,  M.B.  Allen,  H.  Kojouharov,  and  B.  Chen,  ’’Finite-element  solution  of  reaction- 
diffusion  equations  with  advection,”  invited  paper  in  Proceedings,  Computational  Methods 
in  Water  Resources  XI,  ed.  by  A.  Aldama  et  al.,  Cancun,  Mexico,  22-26  July,  1996,  Com¬ 
putational  Mechanics  Publications,  Southampton,  U.K. 
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We  accomplished  the  following  tasks: 


We  developed  a  one-dimensional  finite-element  model  of  solute  transport,  including 
the  effects  of  equilibrium  adsorption  via  a  user-prescribed  isotherm.  The  isotherm 
parameters  can  vary  spatially. 

We  developed  a  methodology  for  assigning  random  isotherm  parameters  (i.e.  the  ad¬ 
sorption  field)  according  to  a  user-prescribed  probability  distribution. 

We  made  a  large  number  of  runs,  based  on  different  statistical  realizations  of  the 
random  adsorption  fields,  to  test  the  hypothesis  that  random  adsorption  has  a  macro- 
scopically  dispersive  effect,  at  least  when  the  isotherm  is  linear.  (Current  work  focuses 
on  quantifying  the  apparent  dispersion  and  extending  the  results  to  more  realistic, 

nonlinear  isotherms.) 

We  began  writing  a  modular  code  to  solve  the  transport-adsorption  equations  m  two 
space  dimensions,  where  the  performance  of  a  numerical  code  is  much  more  sensitive 
to  the  choice  of  numerical  method  used. 


